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We observed the expansion of vortex-free Bose-condensates 
after their sudden release from a slowly rotating anisotropic 
trap. Our results show clear experimental evidence of the 
irrotational flow expected for a superfluid. The expansion 
from a rotating trap has strong features associated with the 
superfluid nature of a Bose-condensate, namely that the con- 
densate cannot at any point be cylindrically symmetric with 
respect to the axis of rotation since such a wavefunction can- 
not possess angular momentum. Consequently, an initially 
rotating condensate expands in a distinctively different way 
to one released from a static trap. We report measurements 
of this phenomenon in absorption images of the condensate 
taken along the direction of the rotation axis. 

PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj, 67.90.-Hz 



In their recent theoretical paper, Edwards et al. cal- 
culated the time evolution of the expansion of a slowly ro- 
tating BEC. At higher rotation rates, greater than about 
half the radial trap frequency, the superfluid nature of the 
condensate has been demonstrated strikingly in the for- 
mation of quantised vortices ||^ . However they point out 
that the superfluidity of the condensate, together with 
the conservation of angular moment, also leads to clearly 
observable consequences at lower rotation rates, under 
conditions where no vortices are present. We report here 
the experimental veriflcation of these predictions by time- 
of-flight measurements on a condensate released from a 
rotating anisotropic trap ||^. These data give further 
evidence for the superfluidity of Bose-condensed gases. 
Other recent experiments on this important aspect of 
Bose-condensation include the observation of quantized 
vortices [ 1 ,D , the absence of dissipation below a critical 
velocity [5| and the quenching of the moment of inertia 
in the scissors mode |6|. 

The solutions of the hydrodynamic equations for su- 
perfluids (that also describe the scissors mode 0|) show 
that a condensate always has a moment inertia less than 
that of a rigid body of the same mass distribution 
For the special case of cylindrical symmetry, the super- 
fluid has zero moment of inertia about the symmetry 
axis. Thus a condensate with some angular momentum 
can never have a circular shape when viewed along the 
rotation axis, since this would imply the unphysical sit- 
uation of infinite angular velocity. An initially rotating 
condensate therefore expands in a different way to one 



released from static trap. We observed this phenomenon 
in absorption images of the condensate taken along the 
direction of the rotational axis - the projection of the 
cloud in this direction never becomes circular (an aspect 
ratio of unity). The evolution of the condensate density 
distribution is calculated using the same hydrodynamic 
equations as in 0, and agrees well with our data. We 
give a brief review of the underlying theoretical aspects 
before describing the experimental procedure and results. 

In our experiment we have an anisotropic harmonic 
potential with three angular frequencies uJx < ujy < ojz- 
The potential rotates about the z-axis with angular fre- 
quency r2. It is convenient to define two new parameters: 
the trap deformation, A (w^ — L0y)/{uj1 + ujy) and the 
mean of the frequencies in the plane of rotation defined 
as = (w^ +LUy)/2. In the hydrodynamic limit, a con- 
densate rotating in this trap displays a quadrupolar fiow 
pattern, as described in fig.0, and has a wavefunction 
which looks like: 



(1) 



where p{r) is the condensate density given in EqJ^ be- 
low. The quadrupolc frequency, v, can be found from 
the solutions of the following cubic equation M: 
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where we introduced the dimensionless quantities D = 

V juJm and f2 = ^jiOra- 

In a rotating trap, the condensate experiences effec- 
tively different trapping frequencies due to its quadrupo- 
lar motion and the condensate density can be written 
as: 



p(r) = - \-—(uj^x +ujyy ) 
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The modified frequencies uJx and ujy are given by 
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where f is found from Eq.g. 

Hence, the aspect ratio of the condensate in the trap 
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where Rx , Ry are the sizes of the condensate in the x and 
y-direction, respectively. The effective chemical potential 
and condensate sizes in the trap can also be calculated 
from EqsJ^. 

To calculate the expansion of a condensate when re- 
leased from the anisotropic harmonic potential, we use 
the following ansatz of spheroidal form for the conden- 
sate density p{r,t) and the velocity field, v(r,i), (as in 

0): 

p(r, t) ^ po{t) - px{t)x^ - Py{t)y^ - Pz{t)z'^ ~ Pxy{t)xy (6) 
v(r, t) = -V {vx{t)x^ + Vy{t)y^ + v,{t)z'^ + Vxy{t)xy) ij) 

Inserting this ansatz into the hydrodynamic equations of 
superfluids yields a set of nine coupled differential equa- 
tions for the expansion parameters, which we integrate 
numerically. The initial conditions for the density are 
found from Eq.||. Note that the cross-term Pxy(S^) is as- 
sumed to be zero at the instant when the condensate is 
released from the trap. During the evolution in time of 
flight, this term grows and can become comparable in 
size to Px, so that the elliptical condensate rotates with 
respect to its initial release angle. All components of the 
velocity field are initially zero, except for the cross term, 
which is the quadrupole velocity field of the rotating con- 
densate, given by Vxy = 2j/ (see Eq.|l|). The behavior 
of the condensate is thus completely determined by the 
above initial conditions and the nine coupled differential 
equations 0]. 

The angle of the condensate and its sizes can be found 
by diagonalising the quadratic form defined in Eq.|[ One 
finds, in time of flight, that at first the condensate ex- 
pands along its smaller axes until the aspect ratio has 
reached a critical value close to spherical. The condensate 
has a quadrupolar flow pattern and an associated angu- 
lar momentum. Conservation of this angular momentum 
does not allow the condensate to become spherical, as a 
spherical quadrupolar flow pattern has no angular mo- 
mentum. Thus, when approaching this critical value in 
the aspect ratio, the condensate rotates quickly and con- 
tinues to expand along its original long axis. The rapid 
increase in angular velocity close to the critical aspect 
ratio arises from the decrease of the moment of inertia 
when the aspect ratio approaches spherical [Q. However, 
energy conservation limits how fast the condensate can 
rotate and thus the aspect ratio cannot become smaller 
than a certain value. This behaviour is in contrast to that 
of a thermal gas which can become spherical and preserve 
angular momentum, since there are no constraints of ir- 
rotationality. Thus it is essentially the superfluid nature 
of the Bose-condcnsatc which prevents it from becoming 
spherical. 

To observe this behaviour we evaporatively cooled 
^^Rb atoms in a TOP trap to a temperature of 0.5 T/Tc, 
where Tc is the critical temperature for condensation. 



This produced Bose-condensates of 1.5 x 10^ atoms, in 
a trapping potential with frequencies Ux/'ii^ = LOyj^Ti = 
124 Hz and lOz/2tt = 350 Hz. We then made the trap 
elliptical with ujyjujx = 1-4 (corresponding to A = 0.32) 
by changing the ratio of the two TOP-field components 
to BxjBy — 4.2. The eccentric trap was rotated by mod- 
ulating the high frequency sinusoidal TOP signal (7 kHz 
frequency) at the low trap rotation frequency 57, as de- 
scribed in ||]. The condensate was spun up in 500 ms 
during which time the eccentricity was ramped up from 
zero to its final value. Then the condensate was left in the 
rotating trap for another 500 ms before it was released 
at a well defined point in the trap rotation, from which 
the angles of the expanding cloud were measured. Fig- 
ure ^ shows typical absorption images of the condensate, 
taken along the axis of rotation, after different expansion 
times. For these pictures the trap was rotating at 28 Hz 
and at the instant of release the long axis of the cloud 
was horizontal. The angle and aspect ratio of the cloud 
was obtained from a 2D parabolic fit to the density dis- 
tribution. Figure ^jb shows the minimum aspect ratio of 
1.31, reached after 4 ms of expansion. The angle in- 
creases steadily to an asymptotic value of '-^ 55 degrees 
in fig.|d. 

It was necessary to go to a high trap eccentricity so 
that a clear deformation of the cloud could be observed. 
We investigated the evolution of the condensate in time 
of flight at two different trap rotation frequencies r2/27r = 
20 and 28 Hz. In both cases the condensate was released 
from a trap of l^JxI'^t^ — 60 Hz, uiy/2'K = 1.4 x 60 Hz and 
ijJz/2tt — 206 Hz. Fig.|| shows the calculated change of 
angle of a condensate in time of ffight after being released 
from a trap rotating at f2/27r = 20 Hz (solid line) and 
Q,/2tt — 28 Hz (dotted line), with the experimental data 
points superimposed. In both cases the angle evolves in 
a similar manner and reaches 45 degrees after about 6 
ms. After 18 ms the angle assumes an asymptotic value 
between 55 and 60 degrees. 

Figure ^ shows the theoretical prediction for the evo- 
lution of the condensate's aspect ratio in time of flight, 
with experimentally measured values superimposed, for 
0/27r = 28 Hz (upper curve), for f7/27r = 20 Hz (mid- 
dle curve) and for release from a nonrotating trap (lower 
curve). The data clearly demonstrate how the aspect ra- 
tio for an initially rotating condensate decreases up to a 
critical point, which is reached after approximately 4 ms. 
From that point on it does not continue to expand along 
its minor axis but the aspect ratio increases again because 
the condensate cannot become circular under these con- 
ditions. However, the condensate released from a static 
trap has no velocity field which prevents it from becom- 
ing circular and the aspect ratio decreases steadily from 
Ljy/ojx to a final value of less than one. Every experi- 
mental point displayed is the average of several measure- 
ments. For each time of flight we had to refocus our 
imaging system as the atoms move out of focus under 
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gravity when released from the trap. However incorrect 
focusing would only result in a more circular image and 
the minimum value for the measurement of the aspect 
ratio of rotating condensates was never consistent with 
unity. There is remarkable agreement between the ex- 
perimental data and the theoretical predictions. How- 
ever, we observed a small deviation of the experimental 
data from the predicted values for the higher rotation 
frequency (r2/27r — 28 Hz, shown in the upper curve). 
We do not know the origin of that deviation. 

Our results show that an expanding vortex-free Bose- 
condensate with some angular momentum refuses to be- 
come circular about the axis of rotation, as predicted 
by Edwards et al 0. This provides direct evidence 
that Bose-condensed gases flow irrotationally as a conse- 
quence of their superfluidity. These measurements com- 
plement the observations of the scissors mode |^ where 
the irrotational flow was deduced from measurements of 
the oscillation frequencies of a trapped condensate - a su- 
perfluid has less moment of inertia than a classical object, 
of the same size, undergoing rigid body motion hence the 
superfluid oscillates faster. Other measurements of the 
flow of a Bose-condensed gas would be of great interest, 
for example flow through a narrow tube ]To| , or array of 
holes, analogous to the superleak experiments with liquid 
helium. 
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FIG. 1. The velocity field of a vortex-free condensate in 
an elliptical trap, rotating at frequency Q. The small arrows 
of the velocity field indicate the direction of the quadrupolar 
flow. 



[1] M. Edwa rds, C.W. Clark, P. Pedri, L. Pitaevskii, and S. 

Stringari |3ond-mat/010608C 
[2] K.W. Madison, F. Chevy, W. Wohlleben, and J. Dal- 

ibard, Phys. Rev. Lett. 84, 806 (2000); J.R. Abo-Shaeer, 

C. Raman, J.M. Vogels and W. Ketterle, Science Express 

(online), March 2001 
[3] J. Arlt, O.M. Marago, E. Hodby, S. A. Hopkins, G. 

Hechenblaikner, S. Webster and C. J. Foot, J. Phys. B, 

32, 5861 (1999) 
[4] M.R. Matthews, B.P. Anderson, PC. Haljan, D.S. HaU, 

C.E. Wieman, and E.A. Cornell, Phys. Rev. Lett. 83, 

2498 (1999) 

[5] C. Raman, M. Kohl, R. Onofrio, D.S. Durfee, C.E. Kuk- 

lewich, Z. Hadzibabic, and W. Ketterle, Phys. Rev. Lett. 

83, 2502 (1999) 
[6] O.M. Marago, S.A. Hopkins, J. Arlt, E. Hodby, G. 

Hechenblaikner and C.J. Foot, Phys. Rev. Lett. 84, 2056 

(2000) 



3 



(a) 



(b) 



(c) 



(d) 



30|im 







FIG. 3. The angle of the condensate plotted against the 
time of flight. The open circles and the filled squares de- 
note the measured angles after release from a trap rotating 
at Q./2ii — 20 and Q./2n — 28 Hz, respectively. The solid and 
the dotted line are the respective theoretical calculations. 
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FIG. 4. The aspect ratio of a condensate in time of flight. 
Initially rotating condensates (upper and middle theoretical 
curves) exhibit a strong backbending effect after 4 ms and 
the condensate never becomes circular. However, after release 
from the non-rotating trap, the aspect ratio decreases steadily 
and inverts (lower curve); it is unity at about 6 ms. 
Upper curve and filled squares, i}/2n = 28 Hz 
Middle curve and open circles, = 20 Hz 

Lower curve and filled triangles. Static trap. 



FIG. 2. Typical images of the condensate at different times 
after release from a trap rotating at 28Hz - (a) 1.09 ms (b) 
3.97 ms (c) 12.00 ms (d) 16.06 ms. At the instant of release 
the long axis of the cloud was horizontal. 
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